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Abstract 

We examine relation between neutrino oscillation parameters and prediction of lepton 
flavor violation, in light of deviations from tri-bimaximal mixing. Our study shows that 
upcoming experimental searches for lepton flavor violation process can provide useful 
implications for neutrino mass spectrum and mixing angles. With simple structure 
of heavy right-handed neutrino and supersymmetry breaking sectors, the discovery of 
t — > ^7 decay determines neutrino mass hierarchy if large (order 0.1) reactor angle is 
established. 
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1 Introduction 



In the last decades, one of the most striking developments in particle physics beyond the 
standard model (SM) is the experimental establishment [1] of neutrino masses and the large 
mixing property, which is quite different from the small mixing in the quark sector. Neutrino 
oscillation experiments have revealed neutrino mass-squared differences and its flavor mixing 
angles. Notably, a simple form of mixing matrix, referred to as tri-bimaximal mixing, is well 
descriptive of the observed mixing structure [2]. Vast numbers of flavor models have been 
proposed in order to derive the tri-bimaximal mixing [3]; thus, from a theoretical viewpoint, it 
is one of the most important subjects to realize difference between tri-bimaximal and observed 
mixing angles. In addition, recent results of the three- flavor global data analysis [4,5] indicate 
non-zero #13; as the best-fit value, not so small one sin $13 ~ 0(0.1) is obtained. Therefore, 
it seems interesting to examine deviations of neutrino mixing angles from the tri-bimaximal 
pattern [6], which leads to sin 2 12 = 1/3, sin 2 #23 = 1/2, and sin 2 13 = 0, for the coming 
future precise experiments. 

The present knowledge of neutrino parameters (i.e. masses, mixing angles and phases) is 
not only important for low-energy characters of neutrinos, but also a key ingredient of the 
origin of flavor structure in the SM fermions. It is thus meaningful to make clear the relation 
between neutrino parameters and high-energy phenomenologies. Among them, charged lepton 
flavor violation (LFV) process would give an intriguing clue, since supersymmetry (SUSY) 
and the seesaw mechanism [7] could enhance LFV as reachable in near future experiments [8]. 
The seesaw mechanism naturally provides desired neutrino mass scale, and predicted LFV 
fractions are affected by low-energy neutrino parameters and heavy Majorana masses via 
SUSY breaking terms [9, 10]. 

In this Letter, we examine the relation between neutrino parameters and LFV predic- 
tion, in light of the tri-bimaximal mixing and the recent precision oscillation data. We use 
a particular parametrization [11] of the MNS matrix [12], where the tri-bimaximal mixing 
is taken as its zeroth order approximation, and give a detailed analysis of LFV prediction 
especially for a simple case of right-handed Majorana mass and SUSY breaking structures. 
Our study shows that upcoming experimental LFV searches provide useful implications for 
neutrino mass spectrum and mixing angles. If we possess simple structure of heavy Majorana 
masses, future discovery of r — > /ry implies that inverted hierarchy (IH) and quasi-degenerate 
(QD) neutrino mass spectra are inconsistent with a large reactor angle of order 0.1, though 
normal hierarchy (NH) is still allowed. 

The Letter is organized as follows. In Section [2], we introduce a useful parametrization [11] 
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of the MNS matrix for examining how the mixing angles deviate from the tri-bimaximal ones. 
In Sections[3]andHl relation between neutrino parameters and LFV prediction in a literature of 
the minimal supersymmetric standard model (MSSM) with the seesaw mechanism is studied. 
Detailed analysis is given in Section H] by using the parametrization. Section [5] is devoted 
to study implications for neutrino parameters with future LFV searches. We summarize our 
result in Section [6j 

2 Neutrino parameters and tri-bimaximal mixing 

Current data obtained by neutrino oscillation experiments is consistent with a simple mixing 
structure called tri-bimaximal mixing [2] . We adopt a particular parametrization proposed in 
Ref. [11] to describe the MNS matrix by deviations from exact tri-bimaximal mixing angles. 
It is helpful to systematically analyze the neutrino tri-bimaximality. 

In the basis where the charged lepton mass matrix is diagonal, neutrino mass matrix is 
given by M v = UT> m U T , where D m = diag(mi,m 2 ,m 3 ) with positive neutrino mass eigenval- 
ues 77^1,2,3 and U is the unitary lepton mixing matrix. Three mixing angles and phases are 
involved in the matrix U; using the standard parametrization [13], it can be expressed as 

— S 12C23 — c 12'S23 s 13 e * <5 C12C23 — Si2S23 s 13 el<5 S23 c 13 JP/Vf, (2-1) 

S12S23 - c 12 c 23 s 13 e tS -C12S23 - s 12 c 23 s 13 e lS c 23 c 13 j 

where Cjj = cos 6^-, s^- = sin 6^-, and 6 is the Dirac CP violating phase. Pm stands for the 
diagonal phase matrix which involves two Majorana phases. 

Recent progress in neutrino experiments greatly increase data for neutrino masses and 
mixing angles. The updated result of the three-flavor global data analysis [4] indicates the 
following best-fit values with 3a intervals of three (solar, atmospheric, reactor) mixing angles 
and two (solar, atmospheric) mass squared differences: 

sin 2 # 12 = 0.304±g;gg|, sin 2 6 23 = 0.50+°;^, sin 2 13 = 0.010 (< 0.056), (2.2) 
Am 2 ol = ml -mj = (7.65±g;j») X lO" 5 eV 2 , 

Am 2 tm ee \ml-m\\ = (2.40±°;|) x 10" 3 eV 2 . (2.3) 

The Dirac phase 6 has not yet been constrained by the experimental data. 

The flavor structure could be determined by profound principles such as flavor symmetry 
in a high-energy regime. Although the origin of the flavor is unrevealed, the current neutrino 
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mixing angles are known to be consistent with a simple mixing matrix Utb, 

l ( 2 y/2 
Utb = (^tb)2 3 (^tb)i2 = -7= -1 \/2 \/3 I , (2.4) 



1 -V2 v^y 



where two (non-diagonal) rotational matrices are given by 

A \ / 2/V6 1/V3 0\ 

(^tb) 2 3 = 1/V2 1/V2 , (^ TB )i2 = -1/V3 2/^6 . (2.5) 
\0 -1/^2 l/v 7 ^/ V V 

Although the matrix Utb approximately describes the MNS matrix ( 12. ip . it is an interesting 
subject to investigate the difference between them in reality. For this purpose, the following 
unitary matrix is useful to parametrize the MNS mixing structure: 

(c x c z s x c z s z c \ 

-s x c y - c x s y s z e t5 c x c y - s x s y s z e iS s y c z (IZtb ) 12 , (2.6) 
Sx^y C X CyS Z C C X Sy S X CyS Z & ^yC Z J 

where c w and s w denote cose w and sine„, (w = x,y,z), respectively!*] In the limit where 
e x, y ,z — > 0, U in (12.61) goes back to Utb except Pm- With this parametrization, experimental 
data ( 12.21) indicates 

- 0.092 < e x < 0.038, -0.14 < e y < 0.17, \e z \ < 0.239, (2.7) 

in its 3cr ranges. Since the deviation parameters are much suppressed than 0(1), expansions 
around e x ,y,z — could give a good approximation to physical quantities related to lepton 
mixing angles. 

3 LFV in MSSM with type-I seesaw 

Let us study LFV prediction in MSSM with heavy right-handed Majorana neutrinos for the 
type-I seesaw mechanism [7]. The relevant part of the MSSM superpotential is given by 

Wiepton = L i (Y e ) ij e j H d + LitY^ijVjHu + ^ViiM^ijVj, (3.1) 

where Y e and Y v are charged lepton and neutrino Yukawa matrices. Superfields Lj, ej, Vi 
(i = 1,2,3) and H u ^) include lepton doublets, charged leptons, right-handed neutrinos and up 



■"■One can easily find that e z = #13 from (|2 . 1 1) and (|2.6p . owing to #13 = in the tri-bimaximal limit. 
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(down)-type Higgs doublet, respectively. Mr gives Majorana mass matrix of the right-handed 
neutrinos, whose scale is assumed to be much larger than the electroweak scale (~ 10 2 GeV). 
It is also noted that scale of M R is required to be smaller than the grand unified theory (GUT) 
scale in order to reproduce known solar and atmospheric neutrino mass scales unless Y v is 
much larger than 0(1). Without loss of generality, we take a basis where Mr is diagonal. 

SUSY breaking should be incorporated in realistic models, since it is not exact symmetry 
in Nature. We thus introduce SUSY breaking terms, which in general could be new sources 
of the flavor violation. Although several breaking scenarios have been proposed, one of the 
most economical and predictive ansatz is to assume the universal form of soft terms in a 
high-energy regime. In this case, universal SUSY breaking parameters are listed as scalar 
mass mo, trilinear coupling Aq, gaugino mass M1/2 and the Higgs bilinear coupling. We refer 
to these SUSY breaking parameters as their GUT scale values. 

The flavor violation in the supersymmetric sector is transmitted to SUSY breaking sec- 
tor thorough renormalization group (RG) evolution between GUT and heavy Majorana mass 
scales. At a low-energy regime, one-loop diagrams with SUSY particles give leading correc- 
tions to the LFV processes; the branching fractions are approximately written as 



8(^-4 + 7) 



a 3 



G 2 F m% 



3m 2 + Al 
_8TT 2 v 2 H sm 2 (3_ 



2 

2fllD 12 



tan 2 ^|5y-r, (3.2) 



where a and Gf are the fine-structure and the Fermi coupling constants, vh — 174 GeV, and 
f3 parametrizes the ratio between vacuum expectation values of the two Higgs scalars in H u ^- 
The mass parameter m§ is a typical mass scale of SUSY particles. Note that flavor indices 
only appear in as 

^ Jl/T 

B l3 = v %sin 2 py2(Y v ) ik (Yj) kj ln-^, (3.3) 
t^i MRk 

where M Rk denotes the fc-th eigenvalue of Mr, and Mq ~ 10 16 GeV is the GUT scale. Thus 
the flavor dependence in LFV branching fractions is completely involved in (13.31) . that is 
SUSY breaking parameters do not affect the flavor structure in the approximate formula of 
LFV prediction (ET2j) . 

Although it is generally difficult to reconstruct the combination of neutrino Yukawa matrix 
in (13.31) from low-energy data, if heavy Majorana neutrinos have an approximately degenerate 
mass Mjj and there are no large CP phases except the Dirac phase 5 in the neutrino sector, 
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the branching fractions (13. 2\i are tightly connected to neutrino parameters as 



„3 



ii + l) 



n- ?>m 2 + A l 



ln^) |M», (3.4) 

3 

% = ^™ fc (tn^ T w (3.5) 

k=l 

Note that Pm in (12.61) disappears in the factor 6y. Here the low-energy neutrino parameters, 
namely their masses, mixing angles and the Dirac phase, are completely involved in (I3.5P ; thus 
the LFV branching fractions depend on SUSY parameters in a flavor independent manner. 
We will discuss more general cases with non-degenerate Majorana masses in the end of the 
next sect ion 

4 LFV prediction around tri-bimaximal mixing 

In this section, we discuss the relation between neutrino parameters and prediction of LFV 
processes. The parametrization (12.61) allows us to handle deviations from tri-bimaximal mixing 
in a systematic way. In a definite framework for high-energy theory, we analyze how LFV 
prediction depends on small deviation parameters with current neutrino oscillation data. 

Analytical results Applying the parametrization (12. 61) to describe the LFV prediction, b{j 
is explicitly written as 

it) 

m i2 r e 



){2\[2c 2x + s 2x ) + ^^c z s z (c y + s y ) |3m 123 + m 12 (c 2x - 2^2s 2x ) 



613 = --^^c z (c y + s y )(2V2c 2x + s 2x ) + -^-j=c z s z (c y - s y ) 3m 123 + m 12 ( C2 X - 2 V2s 2x ) 

hs = ^cl(cl-sl)-^(l + sl)(cl-sl)(c 2x -2V2s 2x ) 

~—m 12 s z [{c y - s y ) 2 - (cy + s y ) 2 e~ 2iS ] (2\[2c 2x + s 2x ), (4.1) 

where c 2x = cos2e :r and s 2x = sin2e x . Note that neutrino masses appear just as particular 
combinations 

777,12 = m 2 -m 1 , m 123 = 2m z -m 2 -m x . (4.2) 



§ One can also discuss relation between neutrino parameters and LFV prediction with other seesaw mech- 
anisms than the conventional type-I [14]. For instance, in the type-II seesaw scenario [15], neutrino masses in 
bij (|3.5p are replaced by those squared, and mixing parameter dependence of LFV prediction can be studied 
as our analysis. 
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It is also noted that bi 2 (e y , e z ; 5) = —b 13 (—e y , —e z ; 5) = —bi 3 (—e y , e z ; 5 ± 7r) holds between 6 12 
and fe 13 ; the relation is just the /i-r symmetric property in the tri-bimaximal limit. 

As mentioned in Section [21 the observed values of lepton mixing angles are consistent with 
the tri-bimaximal mixing pattern. Focusing on the LFV prediction with the tri-bimaximal 
limit (e x ,e y ,e z ) — > (0,0,0) in (14.11) . we obtain 



m 12 

9 ' 

2 ^12 



B(n^ei) oc I&12I 2 -> (-^) 

B{r^ ei ) oc |6 13 | 2 - ^fj' 

K( \ 17 1 2 1 / ^i2\ 2 m 2 23 mi 23 mi2 m 2 2 

B(r^ 7 ) oc |6 23 | --^-—j =— _ + _, (4.3) 



m\ 



where only | &23 1 2 depends on mi 23 and I&12I 2 = I&13I 2 holds. 

Experimental values of solar and atmospheric neutrino mass differences indicate that mi2 3 
is much larger than m 12 . The ratio m = m^/m^ depends on the neutrino mass spectrum. 
For NH of neutrino masses (mi < m 2 < m s ), we obtain 



m = ^ AmL + f~ m ^ > - J^-, (4-4) 

2v/Am 2 tm + m\ - ^ Am 2 soX + m\-m x 2 \\ Am 2 t 



while for the IH case (m 3 < mi < m 2 ), we have 



^ = VAml tm + Am s 2 ol + mf - ^ Amj tm + m\ _^ 1 Am 2 ol ( j _ 

2m 3 - VAmL + Am 2 ol + m 2 - ^/Am^ + ™§ 4Am 2 tm ' 



The last expressions in (I4.4p and (14.51) imply the massless limit of the lightest neutrino, where 
Am 2 ol /Am 2 tm ~ 0.04 with the present data. We plot in Fig. [I] the ratio m as the function of 
the lightest neutrino mass eigenvalue m re f. It shows that the ratio becomes relatively large 
(small) if the neutrino mass spectrum is NH (IH). With being large value of m re f, namely QD 
mass spectrum limit, m takes similar values for normal and inverse ordering cases. 

In the tri-bimaximal limit, B([i — > cy) and B(r — > ey) in H4.3[) are much suppressed than 
Bir — > fi'y) since they do not involve mi2 3 . This implies that these processes are sensitive to 
deviations from the tri-bimaximal mixing. As argued in Section [2j the deviation parameters 
in (12. 6p are sufficiently small, so that we can use them as expansion parameters in LFV 
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Figure 1: Absolute values of the ratio between m 12 and m 123 are shown as functions of 
reference neutrino mass scale m re f; for the NH (IH) case we take m re f = mi(m 3 ). Colored 
bands indicate the predicted values with 3a ranges of input parameters Am 2 ol and Am 2 tm . 
Lines in the bands correspond to the plot with the best-fit values in (12. 3p . 



branching fractions. Up to 0(e 2 ), one can obtain the following expressions: 

|6 12 | 2 ~ fn\ 2 + V2mj 2 e x - 2m\ 2 e y + ^j^(m 12 + mi 23 )e 2 cos<5 - ^-fh\ 2 t 2 x - 2\^m 2 l2 e x e y 

V 2 2 

-^( 7r "i2 - 2mi 2 mi23 - m 2 l23 )e 2 z - ^-(7m 12 - m 123 )e x e z cos 5, (4.6) 
o 2 

I t 1 2 1 / ~ \ 2 ^^12 / ~ \ ^^12 (ft ~ \ 9 1 / ~ \22 

I &23 1 - Y^(^i2-m 123 ) - — (mi2 -mi 23 )ex + — (7mi2 + 771123)^ - -(m i2 -m 123 ) e y 

+ 0(^12 _ m i23 + 16m^ 2 sin 2 5)e 2 z - V2m 12 (m 12 - m 123 )e y e z cos 5, (4.7) 
o 

where rh\ 2 =mi 2 /3 and | &13 1 2 (e^, e z cos 5) = | £> 12 1 2 ( — e y , — e z cos 5). 

From the fact that m\ 23 is much larger than m\ 2 , the typical correlation between neutrino 
parameters and the branching fractions can be understood. In (14. 6p . the terms which remain 
in the tri-bimaximal limit are proportional to m 2 2 , and mi23 always appears with involving 
the deviation parameters, especially e z . Thus | &i2(i3) 1 2 , namely £>(/i(r) — > ej), is sensitive to 
the parameters, while B(r — > yrf) does not have so large dependence on them since m 2 23 is 
the dominant term in \b 23 \ 2 . Moreover, the leading contributions in I&12I 2 and 1 623 1 2 can be 
expressed as 

M 2 - ^(™ 2 + ^rne z cos5 + 9 -el ) )+... } \b 23 \ 2 ~ % + •••, (4.8) 
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since the ratio m is also the small quantity as well as deviation parameters. One can find that 
e z plays an important role to determine B(fi(r) — > ej) and that e x and e y are less effective 
because they only appear as sub-leading corrections. The e z dependence is controlled by the 
neutrino mass spectrum and the Dirac phase through m and cos 5. For example, I&12I 2 is 
minimized at 

2\/2 

e z = #13 ~ mcos<5 ~ — mcos<5. (4.9) 

3 

Note that the value of e z in (14. 9p has significant distinction between NH and IH by the 
magnitude of m. We can check these properties with numerical analysis. 

Numerical searches We discussed the neutrino parameter dependence of the LFV predic- 
tion using the approximate formula in (13. 2\\ . To make our study more complete, we proceed to 
numerical examination of the LFV prediction. Here we take SPSl a [16] for SUSY particle mass 
spectrum; SUSY breaking parameters are fixed as (m , M1/2, A ) — (100,250, —100) GeV at 
the GUT scale and tan (3 = 10. SUSY parameter dependence is mostly flavor blind, and has 
been greatly studied [8]. If one takes other types of SUSY mass spectrum, following results 
do not alter as long as the universality of the SUSY breaking is assumed. 

Given set of SUSY parameters, we numerically estimate RG evolutions between GUT and 
electroweak scale taking heavy right-handed neutrinos into account. Above the right-handed 
neutrino mass scale, which is taken as Mjj = 10 14 GeV in the analysis, the right-handed 
neutrinos are decoupled with the theory. Two-loop RG equations for gauge and Yukawa 
couplings, and one-loop ones for the soft SUSY breaking parameters are numerically solved. 
LFV fractions are estimated with one-loop diagrams in the SUSY particle mass eigenbasis 
rather than the mass-insertion approximation. 

As stressed in the analytical discussion, among the three deviation parameters in (12. 7p . e z 
is crucial for the prediction of n(r) — > ey branching fractions. To see its dependence, we plot 
the LFV predictions as the functions of e z in Fig.[2j The prediction of B(t — > 717) is insensitive 
to size of e z , by contrast that of B(fi(r) — > ej) strongly depends on. The branching fraction 
B(fi — > ej) is minimized around e z ~ — 0.1(+0.01) for the case with NH (IH and QD) mass 
spectrum, as shown in (I4.9p . These results are consistent with the previous argument using the 
analytical expressions. Note that £>(/i(r) — > ej) is highly suppressed than B{r — > 717) around 
the minima. This is an important point to extract implications for neutrino parameters from 
future LFV searches, and we discuss the issue in the next section. 
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Figure 2: Predictions of LFV branching fractions as the functions of e z are shown. Black, 
red and gray plots correspond to £>(/i — > ej), B(r — > fij) and £>(r — > ej), respectively. 
Each figure corresponds to different neutrino mass spectrum: from left to right figures, NH 
(m 1 = 1CT 4 eV), IH (m 3 = 10~ 4 eV), and QD (mi = 10 _1 eV) mass spectra are taken, 
respectively. The mass squared differences are fixed to central values in (12.31) and the Dirac 
phase is taken as 5 = 0. e x and e y are scanned with 3a ranges in (12. 7p . 



Effects of heavy Majorana mass hierarchies In general, heavy Majorana masses have 
non-degeneracy and it affects LFV prediction. In order to incorporate the non-degeneracy of 
Mr, it is convenient to focus on in (13.31) rather than by in (13.51) . By can be generally 
rewritten as follows [9]: 

3 M 
Bio = Y.( U * V ^ RV ^^ V ^ R]V V^ uT )^J^i ( 4 - 10 ) 
t^i MRk 

where = diag{^m{, y/rn^, yjm£), V^j^ = diag(y/M R1 , \JM R2 , V M R3) and a complex 

matrix R satisfies RR T = 1. The additional mixing matrix R appears in Y V Y^ because 
the right-handed neutrino mixing is not unphysical. As a result, LFV prediction generally 
depends on the mixing structure of R, as minutely studied in [9]. 

If the right-handed mixing is not important for low-energy neutrino parameters, namely 
R ~ 1, B^ is simplified. Especially R — > 1 leads to 

B tj ~ J2M k (U*) tk (U T ) kj , Mi = mi M m \n— (4.11) 

It is obvious that B^ has a similar form to b^, where low-energy neutrino mass rrii in (13. 5p is 
replaced by a combination of light and heavy neutrino masses. Hence, explicit expression of 
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Bij is easily obtained by (14.11) . In this case, LFV prediction depends on both neutrino masses 
through A^i2 = A^2 _ -Mi and M.%23 = — M. 2 — -M-i- in particular, the new mass ratio 
M = M 12 /M 123 is essential to determine the neutrino parameter dependence. With typical 
mass hierarchies of Mr, LFV prediction has examined in Ref. [17]. Corrections to R — 1 also 
affect LFV prediction when right-handed neutrinos have large non-degeneracy. For example, 
if R matrix in (I4.10p contains a small mixing angle ^0 then \i — > ey prediction is modified 
with including the counterpart of (14.81) ; up to O(Kij) contribution except for 0(Kij ■ € x ,y,z) 
terms, the prediction is explicitly written as follows: 



\B 



12 



2 _ yVt 123 



M 2 ^ L-.2 , 3 - . 9 2 c V2A 12 + VQA 23 + 2V3A 13 
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where we use the notation 



M 2 + ^A^e z cos 5+-e 2 z + M- \ A + • 

y/2 8 M l23 



(4.12) 



R = I -«i2 1 «23 J , = «ij y/m~m~ I M Ri log — M ra log ^— 1 . (4.13) 



The contribution due to R ^ 1 strongly depends on heavy Majorana mass hierarchy. One 
can similarly see that the leading contribution from Kij appears in for the r — > ey and 
r — > predictions. 

The parameter dependence of the branching fractions is modified from the degenerate 
heavy Majorana case by mainly given difference between m and M.. Nevertheless, a particular 
Majorana mass spectrum is taken such as is dominantly heavy, then similar discussion 
to the degenerate case is possible as long as effects from R in (I4.10p is sufficiently smallf]] 
Further study on effects of R and Majorana masses from our viewpoint is also important and 
future task. 



5 Probing neutrino parameters with LFV searches 

Finally, we investigate possible implications for neutrino parameters from future LFV searches 
with the analysis obtained in the previous section. In the following we concentrate on a limited 
scenario where R = 1 is assumed, and show how future LFV searches give constraints for 
neutrino mass spectrum and 6*i 3 = e z . 

^Here Ky is assumed to be real. Complex phases of R can bring further modification into LFV prediction [9]. 

"in a class of models where all the leptonic flavor violation originates in the charged lepton sector correspond 
to R = 1 as in Ref. [9]. This is equivalent to the case that neutrino Yukawa and heavy Majorana mass matrices 
can be taken as simultaneously diagonal. It is notified that such the situation is approximately realized in 
some Eq and SO(lQ) grand unified models, called lopsided mass structure [18]. 
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Experimental discovery of lepton rare decay processes £j — > £j + 7 is one of smoking gun 
signals of physics beyond the SM; thus several experiments have been developed to detect 
LFV processes. The present experimental upper bounds are given at 90% CL. [19,20]: 

B(n -> erf) < 1.2 x 10" 11 , B{r -> /j 7 ) < 4.5 x 10~ 8 , B(r -> e 7 ) < 1.2 x 10~ 7 . 

These bounds are to be modified in near future searches. MEG experiment searches fi — > e 7 ; 
the bound is expected to reach £>(/i — > e 7 ) < C(10 -13 ~ 10~ 14 ) [21]. Future B-factories would 
also greatly reduce the r decay upper bounds [22]. In our analysis, we conservatively adopt 

B(fi -> e 7 ) < 10- 13 , B(r - a*7) < 10^ 9 , B(t - e 7 ) < 10" 9 , (5.1) 

as upcoming upper bounds of LFV fractions. Since the bound for B(fi — > e 7 ) is most severe, 
— > e 7 )/jB(r — > £*(e) 7 ) must be sufficiently suppressed as 10~ 2 ~ 10 -5 ^ in order to observe 
both /x and r decay processes. Prediction of i3(r — > e 7 ) is always suppressed than B{j — > /i 7 ) 
in Fig. [21 and we focus on r — > /j 7 between the tau decay processes in the following. 

The LFV fractions depend on Mr and SUSY breaking parameters in a flavor blind way, 
and on the neutrino parameters in a flavor dependent manner. As seen in the previous section, 
predictions of /x — > e 7 and r — > e 7 highly depend on e*, but r — >■ /x 7 is nearly independent of 
it. Thus B(t — > /i 7 ) is mostly determined by Mr and SUSY breaking parameters; namely, the 
universal dependence in LFV fractions can be read from B(t — > /i 7 ). Hence we can express 
B(fi — > e 7 ) using e z and S(r — > /i 7 ). Fig. [3] shows contour plots of B(fi — > e 7 ) as the function 
of e 2 and B(r — > /x 7 ) for the cases with NH and IH. In the figure, shaded regions have already 
been excluded by the current experimental bound for /j — > e 7 , and the current and expected 
bounds for r — > /i 7 are shown by the solid and dotted lines, respectively. 

From the figure, one can realize that future LFV searches give us implications for neutrino 
parameters. For instance, if near future experiments discover both £>(/i — > e 7 ) and B(r — > /i 7 ), 
then e z and neutrino mass spectrum are strongly constrained. In the scenario, on the one hand 
for NH case 1 6*13 1 is close to 0.1, on the other hands for IH and QD cases such a large value of 
#13 is not allowed. It is interesting that the above value of reactor angle for NH is in accord 
with the best-fit value reported by recent data analyses. Hence, the LFV discovery excludes 
IH and QD mass spectra when the large #13 is established in experiments like T2K [23] and 
Double Chooz [24]. Though in the analysis the Dirac phase is taken as zero, the result is 
preserved if non-zero value of 5 is incorporated. This is because IH and QD spectra are still 
inconsistent with the LFV discovery and the large #13 value, as easily realized from (14. 9p . 

As another scenario, when future LFV searches discover only B{r — > /i 7 ), the above 
discussion is still valid. Fig. H] shows the prediction of B(r — > /i 7 ) as the functions of e z for 
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Figure 3: Contour plots of £>(/x — > ej) as the functions of e z and B(r — > fij) for the cases 
with NH and IH. The other deviation parameters are set to zero, and the neutrino masses 
are fixed to their central values in (I2.3p . The Dirac phase is taken as 5 = 0. Shaded regions 
have already been excluded by the current experimental bound for /x — > ej, and the current 
(expected) bounds for r — > /ry is shown by the solid (dotted) lines. 

the cases with NH and IH. All the plotted points satisfy B(fi — ► erf) < 1.2 x 10~ n . One can 
see that NH and IH require different values of e z to predict B{r — > 7x7) in future discovery 
range. However, constraints on e z and neutrino mass spectrum are weakened, if B(r — > /X7) 
is sufficiently suppressed than the future experimental limit. 

6 Conclusion 

In this Letter, we have examined relation between neutrino parameters and LFV predictions, 
in light of the tri-bimaximal mixing and the recent precision data. By using a particular 
parametrization for the lepton mixing matrix, which is useful to study difference between the 
MNS and tri-bimaximal mixing matrices, we have explicitly showed that the flavor dependence 
in LFV predictions is controlled by deviation parameters and neutrino mass differences. 

In the setup with universal heavy Majorana masses and soft SUSY breaking parameters, 
we have found that e z and the neutrino mass spectrum are important for predictions of 
£>(/x — > e-y) and £>(r — > cy), while e x and e y are less effective to determine the LFV predictions. 
The branching fraction B[j — > /17) is also shown to be insensitive to neutrino parameters. 



12 




-7 



excluded 



future 
sensitivity 



■ too 
-0.2 



-0.1 







0.1 



0.2 



Figure 4: Prediction of Bir — > 7x7) as the functions of for the cases with NH (black circles) 
and IH (gray triangles). All the plotted points satisfy B(fi — > ej) < 1.2 x 10~ n . The deviation 
parameters e x and e y are scanned with 3a ranges (12.71) . and the Dirac phase is taken as 5 = 0. 
The current (expected) bounds for r — > pq is shown by the solid (dotted) lines. 

In addition, we have discussed the effects from heavy Majorana mass structure, namely non- 
degeneracy in right-handed neutrinos and small mixing angles in R matrix. 

We have examined and extracted the possible implications for neutrino parameters from 
upcoming LFV searches. Future discovery of LFV process can give strong constraints on # 13 
with respect to the type of neutrino mass hierarchy as long as effects from R in (14.101) is 
sufficiently small. In particular, r — > /J7 discovery excludes IH and QD neutrino mass spectra 
if I #13 1 ~ 0.1 would be established. 

In general, inclusive studies of the precision neutrino parameters and LFV could give us 
implications for unrevealed issues in the lepton flavor structure, such as the neutrino tri- 
bimaximality. Further investigations of LFV prediction focusing on effects from the phases 
and the Majorana mass structure are our future works. 
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